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Evaluation model of the grey fuzzy on eco-environment vulnerability 
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Abstract: The basic theory and evaluation index system of eco-environment vulnerability were reviewed. Based on the grey theory and 
fuzzy mathematics, a new comprehensive evaluation method from qualitative to quantitative, called grey-fuzzy evaluation, was proposed for 
evaluating eco-environment vulnerability. It was integrated of Association for Healthcare Philanthropy (AHP), grey correlation analysis, 
grey statistics and fuzzy judgment. The constitutional principle and method of the new evaluation method were given and its feasibility and 
effectiveness were proved by the practical example. 
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Introduction 

The ecological environment is defined as the synthesis of the 
climate, landform, soil, hydrology, zoology and botany, human 
activities and so on (Liu 1995). In recent years, the ability of 
coordination ecosystem oneself is descending continuously and 
the survival environment for the human is presenting an increas¬ 
ingly weak trend due to the sharp growth of the population and 
the unreasonable exploitation to the resources. Therefore, the 
research on the vulnerability of the ecological environment has 
become a popular issue in the domain of the resource and envi¬ 
ronment, especially the evaluation of the vulnerability of eco¬ 
logical environment. The object of the research concerning the 
vulnerability of ecological environment is focused on resource 
and environment system that is susceptible to human activities 
and external environment intimidation. The essence of the sys¬ 
tem, which covers the internal structure and external condition 
and also reflects the mutual impact between the system of envi¬ 
ronment resource and the intimidation of the external environ¬ 
ment, is described by the vulnerability of ecological environment. 
The evaluation on vulnerable ecological environment mainly 
involves the evaluation of instability of the ecosystem owing to 
the disproportion of various factors in space and time and the 
impact evaluation of the ecosystem because of human activities 
and the intimidation of the external environment. The evaluation 
of eco-environmental vulnerability is a synthesis method con¬ 
cerning the assessment of population, production, social devel¬ 
opment and the capacity of the whole environment or the degree 
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of capability (Yang 1992). The assessment is the whole evalua¬ 
tion of ecosystem effects that is for the special attributes of the 
various environmental factors and factors combination, which is 
known as the scope of eco-environment vulnerability and the 
evolution tendency. We combined the comprehensive analysis on 
the reason of frail ecosystem and environment receptor, with the 
analysis method of fixed-quantity and semi-quantitative, in order 
to generalize integrant of eco-environment frailty (Li et al. 2002). 
This work not only takes advantage of the quantitative contrast in 
the regions, also provides the basis for the whole treatment to the 
vulnerability of ecological environmental synthesis. 

The basic theory of grey-fuzzy evaluation on 
eco-environment vulnerability 

Considering the complex system with a typical of the incomplete 
information, we use the grey system theory to cope with the ob¬ 
servational data in some layers by taking the non-statistic 
mathematical method (a better practicability under the condition 
less dates and not meeting the request of the system) in order to 
obtain a better understanding of the interior variety of the system 
and correlation mechanisms in a higher level (Deng 1990). 
However, the fuzzy mathematics is one kind of mathematical 
method that is applied to describe and process the fuzzy informa¬ 
tion (Qu et al. 1988). On the other hand, it brings the mathematic 
into an objective fuzzy field and sets up a bridge between the 
traditional classical mathematics and the fuzzy real world. When 
a fuzzy phenomenon is described by taking this method, its es¬ 
sence and rule are revealed. 

The vulnerability of ecological environment has a typical 
fuzziness such as impact factors (the district of the time and 
space on the plant), the weight of assessable indexes and the 
phenomenon and property itself (Liu 1995). In addition, the eco¬ 
system also has the characteristics of the grey system, including 
the imperfect information, the incomplete relative factors of a 
ring, the incomplete-known systematic structure and the incom¬ 
plete understanding of the principle of the system, especially that 
the distribution between value and time of the ecological envi¬ 
ronment clearly exhibits a very fuzzy process. Therefore, on the 
basis of the both characteristics of the ecosystem, this paper 
builds up a grey fuzzy model of the vulnerability of ecosystem 
by making use of the method of fuzzy mathematics combined 
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with grey system theories. 

The evaluation of the grey fuzzy model on 
eco-environment vulnerability 

Grades of vulnerability of ecosystem environment 

Ecosystem environment can be divided into 4 grades according 
to the vulnerable degree of the ecosystem environment. V={the 
most vulnerable v 1, the more vulnerable v 2, the medium v 3, 
light vulnerable v 4}. 

Establishing the evaluation model with gather hiberarchy 

There are control points with the number of “n” in the valuation 
area. Each control point needs to consider the primary factors 
with the number of “m” which includes “K,” (z'=l,2,...,m) be¬ 
longing to the secondary factors. Set up a structure model with 
three layers: target layer (flimsiness), stipulation layer (the first 
class evaluation factor), Index layer (the second class evaluation 
factor). The model structure is shown as Fig. 1. 



Fig. 1 The structure of evaluation model with gather hiberarchy 


education from universities and technical secondary schools and 
the total population, and nutrition indexed by the Engel Coeffi¬ 
cient. 

Quantitative and qualitative indicators are dealt with stan¬ 
dardization by using different methods. 

The excellent process of the quantitative indexes 
The assessable indexes are divided into three types: cost index, 
efficiency index and zone index. The lower the cost index the 
better and the higher the efficiency index the better, such as the 
using rate of resource, satisfaction degree and so on. The zone 
index means a kind of index which takes a certain fixed zone 
belonging to, such as advanced index. According to the classifi¬ 
cation of the indexes, the Qj also is divided into three subclass set 
further, even if: 

e = US Qn0,=0 ij = 1,2, i*j ( 1 ) 

z=l 

where, The Qj (because of the inconsistence of quantity key link 
in the various indexes, the result is not carried on contrast analy¬ 
sis directly, z=l ,2,3) represents respectively cost type set, per¬ 
formance type set and zone index set, and the 0 signifies an 
empty set. 

No method can be used directly for comparative analysis due 
to inconsistent dimensional indicators. Therefore, for different 
types of superscript, we should adopt different non-dimensional 
methods (Yang 1992). 

The transform equation of non-dimension for cost index is as 
follows: 

d max -d.. 

— -— ( 2 ) 

lJ d™ x - dj m 

i= 1,2,..., m;y'=l,2,...,n and OjG Qj 


Evaluation index system 

According to the vulnerable characteristics of ecosystem envi¬ 
ronment, the vulnerable valuation should be a comprehensive 
work for each main factor of ecosystem environment in the dis¬ 
trict. The key of the vulnerable evaluation is to set up a reason¬ 
able system of assessable indexes. 

The evaluation index system can estimate whether the evalua¬ 
tion result is reliable and possible or not. Because of the vulner¬ 
able changes coming from a coefficient difference between natu¬ 
ral factors and artificial factors, the index system should be a 
synthetic system, and the index system chosen should be overall 
contexts and reflect the actual conditions of natural ecological 
environment. 

Natural factors include architectonic, morphologic characteris¬ 
tics, surficial constructional materials, phytocommunity and 
biotic population type; Climatic factors mainly include light, hot, 
water and other factors, which are the material basis for the eco¬ 
logical environment constituents; the conditions and spectrum of 
economic development are presented by the degree of forest 
coverage, annual precipitation, potential evapotranspiration, per 
capita possession of water resources, per capita income, agricul¬ 
tural modernization degree, power farming agricultural acreage 
as well as chief agricultural acreage, the knowledge population 
expressed by the ratio between people who have received higher 


Where, the dj max and the dj mm represent respectively the maxi¬ 
mum and minimum in the allover evaluations 

4,- ,d 2j ,...4 mj }; 4=mir(4,,4 ; ,-4„ ; } ( 3 ) 

The transform equation of non-dimension for efficiency index is 
as follows: 


d„ - d p 

b = —---—— 

11 dp - dj' n 

i = 1 , 2 , ... , m ; j 

QjS Q2 



, 2 , ... , n and 


The transform equation of non-dimension for zone index is: 


h. = 


\ 4 - 4, 

I 1 u 

ma x^-dp,dp-S 2 \ 

d fj < s t 
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z=l,2,...,m;y=l,2,...,n and Q 2 
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where, [S 1? S 2 ] is the best interval for the O index 
Fuzzy quantification qualitative index 

The trapezoidal fuzzy number is used to express the indexes, 
Among these indexes, “bad” corresponds with trapezoidal fuzzy 
number (0, 0, 0, 0, 2), “a little bad” corresponds with the number 
(0, 0, 2, 0, 2, 0, 4), “common” corresponds with the number (0, 3, 
0, 5, 0, 5, 0, 7), “a little good “corresponds with trapezoidal 
fuzzy number (0, 6, 0, 8, 0, 1, 0), “good” corresponds with the 
number (0, 7, 0, 9, 1, 0, 1), and “very good” corresponds with the 
number (0, 8, 1,0, 1,0, 1, 0). For the convenience of calculation, 
the trapezoidal fuzzy numbers are often written as L-R fuzzy 
number, and the relationship is: trapezoidal fuzzy numbers (cp, m, 
n, (3), L-R fuzzy numbers (m, n, y, 5), y = m- cp, 5 = (3 - n, L-R 
fuzzy number M = (a, b, cp, (3), N= (c, d, y, 5). According to the 
literature (Yang 1992), the equation used for the division of 
fuzzy numbers for the trapezoid fuzzy number can be expressed 
as: 


[M IN] « 


a b a8 + dtp by + c/3 
d 9 c’ d(d + 8)’ c(c - y) 


] 



The overall expectations for the trapezoid fuzzy number (A = (cp, 
m, n, (3)) are: 


I 


A 


fa + m + n + j3 f 


V 


A 


J 



Fuzzy grey correlation analysis 

The original time series must be dimensionless. To avoid the 
impact of inconsistent unit to correlation analysis, the starting of 
the various time series should be coincident after treatment, and 

the new time series were respectively recorded as Iq°^ (t) and 
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(9) 


In order to move comparative time sequence toward the refer¬ 
ence sequence, a linear transform should be done, which can 
make the starting and ending value coincide with the reference 
sequence without the change of curve shape. The new sequence 

transformed is called generating sequence, {t\ 


( 0 ) 


( 0-1 


xr(t) = co + \^-\co-0) 

xf \m)- 1 


( 10 ) 


If T represents time, the difference of absolute value between 
x 0 (0) (t) and x 7 (1) (t) can be expressed as A oit =\ x 0 (0) (t) - x 7 (1) (t) |, 
t= l,2,...,m; z'=l,2,...,n. According to the expression, the mini¬ 
mum difference of absolute value can be expressed as 
A = min min X 0 (0) (z) - X (1) (t)\ an d maximum difference of 


absolute value for the moment can be expressed 
as A max = maxmaxX (0) (t) - x/ (f | • Generally, the distinguishing 

it 

coefficient is 0.5. Actually, the correlation coefficient £ oi(t) is 
variable and correlated with X { and other comparative sequences 
X h (h= l,2,...,n). Because the maximum absolute value can reveal 
the whole characteristic of the system, the distinguishing coeffi¬ 
cient D as a maximum should reflect sufficiently the indirect 
impact of various factors to associated space and play a partial 
role in resisting interference. Namely it can weaken the error 
generated from the abnormal observation of the entire associated 
space. 

The method of defining the distinguishing coefficient D is as 
follows (Li et al. 2002; Koczy 1983) 



5m 



ib-xrit) 


(ii) 


If q=A k /A max , the range of 5 is q<5<2 which should meet the 
conditions that A max >3 Ak , q<5<1.5 q; or A max >3 Ak , 1.5 q<5<2 q 
After ascertaining the distinguishing coefficient, the correla¬ 
tion coefficient can be calculated by the Eq. (12), in order to 
focus on the whole correlation space: 

max max | x 0 (t) - x.(t) | ^2) 

x 0 (/ L ) - x. (t) | +8 max max | x 0 (/ L ) - xft) \ 

i t 



The structure of fuzzy judgment matrix 


Based on hierarchical model, the pair factors in the same layer 
are used in a contrast computing involved the above so as to 
constitute a contrast judgment Matrix in the both factors (Zhao 
1986). 

If u 7 (z'=l,2,...,m) is expressed as the first-grade evaluation 
factor, A—C (judgment Matrix) of the first-grade evaluation 
factor is obtained. 


Pc = 


r ll 

r i 2 

V 

' lm 



r 21 
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r 22 
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(13) 
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(i, j= 1,2,... 

,m) represents the relative important 

value 


1V 1/ 

from u 7 to Uj in the first class evaluation factor Pc [; 

Suppose that Ui represents the second-grade evaluation factor, 
C—Ciai (ai= 1,2,...,K 7 ) the judgment Matrix Pc can be obtained 
in the same way. 


Where, (<d, 0) is a range for generating sequence, generally [1, 
x 0 <0) ( m )]. 
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acquired by these steps, and the eigenvector is the weight multi¬ 
tude Ac= (W 1 ,W 2 ,...,Wm) of the first-grade evaluation factor. In 
the same way, the eigenvector of the second-grade evaluation 
(14) factor in judgment matrix can also be acquired. The eigenvector 
is W ci =(W [u W [2 ...,W iKi )T. According to W c =W{Wy (7=1,2,...,m; 
y=l,2,...,A7), we could establish the weight multitude, 
A ci =WfW ci = (W cil ,W ci2 ,...,W ciki ), of the second-grade evaluation 
factor. 


Ascertaining the weight of each element of assessment 


Establishing fuzzy evaluation matrix for each layer 


The eigenvector is calculated by judgment matrix, which is cor¬ 
responding to the maximum eigenvalue and is the weight value 
of relative importance 

Calculate the product per row element M 7 in judgment matrix. 


™_ us 

= n r iM,j = < 15 ) 

7=1 

1 

Calculation of W 7 , it is expressed by ™Jm . , W=M ™ ; 
Normalization of the vector W: 

W=[W 1 ,W 2 ,...,W m ] T , so W=[W!,W 2 ,...,W m ] T is the required 
eigenvector used to calculate the maximum estimate 2 max in 
judgment matrix 


2 


max 


f (PEA _ 1 y ( pw l 
tr nw, W t 



Where, (PW)j is the / element of the vector PW. 



The test of the weight: since the above eigenvector is the re¬ 
quired weight, it is necessary to proceed consistency check for 
the judgment matrix, in order to prove the rationality. 

CR = TT ; CI = 1 1 C “ m ) (18) 

Rl m — 1 


The fuzzy evaluation matrix with the single-factor evaluation can 
be used as an assessment of various grades of vulnerability 
evaluation. 

The index information of every single-factor in the ecological 
environment is limited. As far as the whole process, it is born 
with the characteristic of the grey system. Thus, the rating of 
grey (grey-fuzzy membership) to each single factor in the vari¬ 
ous grades of evaluation is determined by the grey statistical 
analysis of single-factor. The fuzzy relationship between the 
constitutional factor being composed of the grey-fuzzy member¬ 
ship and grades of vulnerability is defined as the evaluation ma¬ 
trix (cr 7 ). The grey-fuzzy membership of eco-environment vul¬ 
nerability scale can be calculated by grey statistical analysis. 

And I , II ,III,...,co: statistical analysis; 1° and 2°,3°...,m°: un¬ 
favorable for statistical indicators; l,2,3,...,n: statistical grey; f b 
f 2 ,...,f n : grey category weighted functions; Real number dy is the 
object (sample) which is the i statistical to the j indicators data. 
dy, i £ { I ,II,...,co}, {l°,2°,...,m°}, D is the matrix with 

elements in dy ( Turksen et al. 1994). 


D = 
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If F is mapping, Opf k (dy) is a calculation of f k (dy), and er 7 is the 
Weights (sequence): grey ash 1, grey ash 2,..., grey ash n. cr=(<7 /7 , 
Oj 2 ,-~, ffjn), {1°.2 }. WhenF: o/ k (dy 0,1], 

K& {1,2,—,n}, i'G{ I ,II,...,co} JS {l°,2°,...,m F(o/ k (d tJ )) is 
the grey statistics that all statistical indicators for the j target. 

. . w (O n CO 

F{o p f k {d ij ))= X f t d„/'Z X Jk , ( 2 °) 

i = l k = 1 i =1 


fk (dy) satisfy: grey <S> e[x 2 »]; 


where, CR is the random consistency ratio of the judgment ma¬ 
trix, CI the common coincidence indicator of the judgment ma¬ 
trix, and RI the average random coincidence indicator of the 
judgment matrix for the judgment matrix 1-9. The value of the 
RI could refer to the literature (Zhao 1986). 

When CR<0.10, the content consistency of judgment matrix is 
reliable, and also indicates that the weight distribution of the 
evaluation factor is reasonable; otherwise, the judgment matrix 
should be adjusted until the content consistency. The eigenvector 
of the first-grade evaluation factor in judgment matrix can be 


/ A U, )-/.,<</, Z ' 7 

y _ -y* 
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f k (d, i ) = L k (d, i ) 

fk(dij) = 1, 

f k {d ij ) = R k {d i] ) 

f k (d y ) = 0, 


d ,j ~ 
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-d, 


V 



d y s[x ,,x 2 ]. 


d.. — x 7 
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dy e[x p x 3 ] 
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[0, X[, x 2 ]; 




A(rf#)=i. 


<4e[0,x, J 


fk( d ij )= R d d ij )=—-^s[xi ,x 2 ] 


(23) 


x 2 — ^ 


tion, the synthetic evaluation is also divided into two levels. The 
first level is the A Ci , (the weight assignment of evaluation factor 
in the second level that is determined by the hierarchy analytic 
process) and the R Ci (the evaluation matrix determined by the 
first level) carries out synthesis operation in order to come out 
the judgment of the level factor, which is called the first level 
fuzzy synthetic evaluation model: 


B ci = A ci • R ci 

"^11 

^12 

^Iki 

= (W cil ,W ci2 ,-,W ciK ). 
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•" &2ki 


_^ki2 

^"ki2 

^kiki. 

= ( B B ••• B 1 

V Cil 5 Ci2 > > JJ CiKi J 





sign: <S>CE[x 2 co]; G^i,x 2 ,x 3 ]; 
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The second level fuzzy synthetic evaluation models: 

AcKW.W,,... ,WJ;B = Ai-Ri = (W,,W 2 „.. ,W m >Ri 

The process is that B b B 2 ,...,B m (the result of the first level 
evaluation in the second level evaluation factors) is arranged in 
longitudinal row to constitute second level evaluation matrix . 

The overall merit model for the second-grade fuzzy is showed 
as follows: 


B = A i *R i =(W „W 2 ,-W J»R, 


= (W „W 2 ,-W J 


B 

J 
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= ( v,,v 2 ,v 3 ,v.) 


(27) 


Evaluation matrix (cr 7 ) is defined as the fuzzy relationship be¬ 
tween the grey-fuzzy membership and vulnerability scale, 
namely constituted fuzzy evaluation matrix of the first-grade 
factor. 
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where Oy (i,j=l,2,...,ki) is the grey-fuzzy membership, which is 
the first evaluation w 7 to the weights of Uj 

Evaluation models on the grey-fuzzy 

The evaluation model of the grey fuzzy is a mathematics tool 
which the synthesis determination is implemented by adopting 
fuzzy mapping and the fuzzy linear substitution. According to the 
numbers of factors involved in the judgment, the fuzzy compre¬ 
hensive evaluation is divided into the single fuzzy synthetic 
evaluation and the multistage fuzzy synthetic evaluation. In addi- 


On the basic principle of the maximum grade in fuzzy mathe¬ 
matics, the vulnerable scale of ecological environment can be 
evaluated to some extent. 

Examples 

Wuchuan County, Inner Mongolia is located in the south of Inner 
Mongolia plateau, north foot of Yinshan Mountain. It has a typi¬ 
cal temperate, semi-arid continental monsoon climate. The an¬ 
nual average temperature is 3.1 °C and the annual average rain¬ 
fall is about 300 mm, mainly concentrating in July to September. 
The main vegetation type is dry grassland, shrub and desert plant, 
and the main agricultural crops are wheat, oat, etc. The economy 
of agriculture and animal husbandry is backward, and desertifi¬ 
cation and soil erosion are increasingly serious in this region. 

Identification of the hierarchical structure and evaluation indexes 

The eco-environment of the staggered areas between agricul¬ 
ture and animal husbandry in the north of China is very complex, 
which the cause involves many factors and is complex. The key 
is to choose evaluation factors. In this paper, 13 indicators were 
collected according to the principle of desirability, concise and 
comparability. 
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The fuzzy matrix 

The fragile indexes of eco-environment at Wuchuan County were 
calculated by the grey correlation, establishing the fuzzy judg¬ 
ment matrix. 

Based on the calculation of the matrix eigenvalue, eigenvalue 
vector and consistency test, the distribution of the total weight 
and subsystems weight can be obtained. A (Total weight)=(0.66, 
0.15, 0.09, 0.10); Al (natural resources subsystem weight)=(0.68, 
0.23, 0.09); A2(environment subsystem weight)=(0.61, 0.25, 
0.08, 0.06); A3 (economy subsystem weight)=(0.77, 0.16, 0.07); 
A4 (social development subsystem weight)=(0.76, 0.12, 0.12). 

Fuzzy evaluation matrix 

Each evaluation grey weight of various factors was obtained by 
analyzing grey single-factor, and grey membership category. 
Evaluation matrix was defined as the fuzzy between the compo¬ 
nent factors and vulnerability grades. The results were shown as 
follows: 

Natural subsystem evaluation matrix 



0.32 0.48 0.12 0.08 
0.26 0.28 0.22 0.24 
0.28 0.40 0.20 0.12 


Respectively bring the evaluation matrix and weight into the 
grey fuzzy evaluation models, we can get: B= (0.54 0.19 0.19 
0.08). According to the principle of largest membership in fuzzy 
mathematics, we can conclude that the ecological environment at 
Wuchuan County, Inner Mongolia belongs to the extremely frag¬ 
ile district. 

Conclusions 

The evaluation model of eco-environment vulnerability which 
uses grey theories and vague mathematics gathers many experts’ 
idea. The model carries forward the advantages of the classic 
method, also controls and makes up its partial weaknesses such 
as people’s mistakes. The practice confirmed that this method 
was effective and reliable for solving vulnerability evaluation of 
eco-environment. Moreover, it also provides an important basis 
for solving the similar evaluation problems, especially involving 
the fuzzy and multi-level grey-fuzzy analysis. 



0.57 0.22 0.18 0.03 
0.43 0.41 0.07 0.07 
0.48 0.16 0.24 0.12 


Environmental subsystem evaluation matrix R 2 



0.68 0.24 0.04 0.04 
0.46 0.32 0.14 0.08 
0.38 0.40 0.16 0.06 
0.39 0.36 0.13 0.12 


Economic subsystem evaluation matrix R 3 



0.46 0.38 0.09 0.07 
0.44 0.32 0.16 0.12 
0.38 0.36 0.10 0.16 
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